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ABSTRACT 

We have identified a new class of Asymptotic Giant Branch (AGB) stars in the Small and Large 
Magellanic Clouds (SMC/LMC) using optical to infrared photometry, lightcurves, and optical spec¬ 
troscopy. The strong dust production and long-period pulsations of these stars indicate that they 
are at the very end of their AGB evolution. Period-mass-radius relations for the fundamental-mode 
pulsators give median current stellar masses of 1.14 M 0 in the LMC and 0.94 M 0 in the SMC (with 
dispersions of 0.21 and 0.18 M 0 , respectively), and models suggest initial masses of <1.5 M 0 and 
<1.25 M 0 , respectively. This new class of stars includes both O-rich and C-rich chemistries, placing 
the limit where dredge-up allows carbon star production below these masses. A high fraction of the 
brightest among them should show S star characteristics indicative of atmospheric C/O ^ 1, and many 
will form O-rich dust prior to their C-rich phase. These stars can be separated from their less-evolved 
counterparts by their characteristically red J — [8] colors. 


1. INTRODUCTION 


1.1. Asymptotic Giant Branch Stars 
Thermally-pulsing Asymptotic Giant Branch (TP- 


near- 

-infrared (NIR) luminosities (|Maraston et al. 

2006 

Melbourne et al. 2012 Melbourne V Boyer 2013 

and 

to the dust budgets of galaxies ( 

e.g., Zhukovska & 

Hen- 

ning 

2013 Schneider et al. 2014; 

McKinnon et al. 2015|. 

The 

Magellanic Clouds are excel 

ent galaxies for investi- 


gating the global effects of TP-AGB stars because they 
are very nearby (<65 kpc), sufficiently massive (0.5- 
1.5x 10 6 MQ tars ), and have high star-formation rates at 
epochs required to form TP-AGB populations in the 
(M 
IB- 


present day (|McConnachie|2012| and references therein) 
oth 


The TP-AGLTstars in both the Large and Small Magel¬ 
lanic Clo uds ( LMC/ SMC) are therefore extensively stud¬ 
ied. In Boyer et al. (20111) (hereafter Bll), we compiled 


a catalog of TP-AGB stars in the MCs using photometry 
from visible, NIR, and IR surveys and noted a distinct 
population of TP-AGB stars with 8-/l m excesses that 
separated from the bulk TP-AGB population in color. 
We dubbed these the anomalous oxygen-rich (aO)-AGB 
stars. In this work, we investigate the nature of these 
stars using additional multi-wavelength and time series 
data. 

Before the onset of thermal pulses, stars in the AGB 
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evolutionary stage are oxygen-rich and their optical to 
NIR spectra are dominated by O-rich molecules such 
as TiO, VO, and H 2 O. Once thermal pulsing begins, 
newly-synthesized elements are dredged-up to the sur¬ 
face. In intermediate-mass TP-AGB stars, the dredged- 
up material includes carbon, which quickly pairs with 
free oxygen to make CO molecules. If the dredge-up is 
efficient enough, the ratio of carbon to oxygen exceeds 
unity and the excess carbon is free to form carbon-rich 
molecule s such as CN and Cg- Stellar evolution mod¬ 
els (e.g., |Karakas et al.||2002| iMarigo et al.||2008| |2013) 
predict a higher dredge-up efficiency as metallicity de- 
creases. This, combined with less free oxygen, results 
in the rapid formation of C stars with fewer dredge-up 
events. Carbon stars are therefore more easily formed in 
metal-poor environments. 

Whether a star remains O-rich or becomes C-rich has 
drastic consequences for its subsequent evolution, per¬ 
haps affecting the chemical evolution of its host galaxy. 
C-rich TP-AGB stars (C-AGB) easily form dust and 
quickly develop strong, radiation-driven winds. O-rich 
TP-AGB stars (O-AGB) can also exhibit large mass-loss 
rates, but in galaxies like the Magellanic Clouds, it is 
the C-rich stars that currently dominate the stellar dust 
production despite having much smaller numbers than 


vasan et al. 


their O-rich counterparts ( Matsuur a et al. 2009 
|2009 Boyer et al. 2012 ~Riebel et al 
B 


Srini- 


2012 ) 


The aO-AGB stars identified in Bll show IR co. 
termediate to the O-AGB and C-AGB population, and 
their nature was therefore unclear. 


1.2. Photometric Classifications of TP-AGB stars 

O- and C-rich molecular absorption features are strong 
at optical and NIR wavelengths, causing the colors of C- 
AGB and O-AGB stars to be drastically different in NIR 
photom etry wit h commonly used filters. Several works 


(e.g., Cioni et al. 2006|) have shown that O-AGB stars 
exhibit J — As ~ 1 mag, extending along a mostly verti¬ 
cal sequence in Figure [l] upwards from the tip of the Red 
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Giant Branch (TRGB). C-rich stars follow a more hori¬ 
zontal sequence in Figure [lj reaching J — K$ ~ 5 mag, 
though very dusty or very cool O-AG B stars (su ch as 


OH/IR stars and metal -rich stars, e.g.; van Loon et al. 
1998 Boyer et al.J|2013 ) can also show these red colors. 


These molecular features therefore make NIR photome¬ 
try an effective, if imperfect, tool for separating C-AGB 
and O-AGB stars when spectra are not available. 

NIR d ata from the Two -Micron All-sky Survey 
(2MASS; |Skrutskie et al.|[2QQ6| for the Small and Large 
Magellanic Clouds (SIVIC/LIVIC) show this clear separa¬ 
tion between the two types of AGB stars in J — K$ color 


e.g. 


Cioni et al. 2006 Blum et al. 


2011). However, Bll found that, w. 


2006 Boyer et al. 
ren combined with 
8 fi m photometry, a third group of stars distinguishes it¬ 
self from both the C-AGB and O-AGB stars by showing 
J — [8] colors that are redder than the O-AGB stars and 
bluer than the C-AGB stars (Fig. T|) . More than 95% of 
these stars are classifie d as O-rich b ased o n the NIR pho¬ 
tometric criteria from[Cioni et al.|| (2006 1) , so we dubbed 
them anomalous O-rich AGB stars (aG-AGB). 
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TABLE 1 

Photometry Summary 


Catalog filter Resolution Refs. 

(jLtm or band) ("/pix) 


2MASS 1.2, 1.6, 2.2 1" 

SAGE 3.6, 4.5, 5.8, 8.0, 24 1722-274 

WISE 3.4, 4.6, 12, 22 671-1270 

AKARI 3.2, 4.1, 7, 11, 15, 24 279-278 

MCPS 0.37 (U),0.44 (B) 172-178 

0.55 (V),0.81 (I) 172-178 


1 

2,3 

4 

5,6 

7 

7 


References. — (1) Skrutskie et al. 

(2006 

1.(2) 

Meixner et al. 

(2006), (3) Gordon et al. 

(2011 

Li£ 

Wright et al. ( 

2010), (5) Ita et al. 12010 

OBJ 

Kato| 

et aT] (2012), (7) Zaritsky et al. (1997). 




their circumstellar envelopes. Second, these stars may be 
a subset of O-AGB stars with particularly strong and/or 
broad silicate emission, enhancing their 8-/im flux. A 
third possibility is that the aO-AGB stars are S-type 
AGB stars, which are transitioning from O-rich to C-rich 
as they evolve. It is impossible to discern whether one 
of these scenarios is correct with only optical to mid- 
IR photometry. Indeed, we cannot even be sure that 
the aO-AGB stars are truly O-rich without additional 
information. 

In this work, we bring together optical spectroscopy, 
optical to mid-IR photometry, and optical pulsation in¬ 
formation to characterize the aO-AGB stars with respect 
to the other O-A GB an d C-AGB stars. One finding, dis¬ 
cussed in Sect ion [2. 3. 1[ is that the aO-AGB stars are not, 
in fact, all O-rich. We therefore do not wish to propagate 
this potentially confusing nomenclature and will instead 
refer to the aO-AGB stars simply as a-AGB stars. While 
the name has changed, we emphasize that the a-AGB 
stars are photometrically classified with a scheme iden¬ 
tical to the aO-AG B classification scheme described by 
Bll (Section |2.1.1 ). 

With the rich data set compiled here, we com pare 
the chemistry (Section pTH , pulsation (Section pL2| ), stel¬ 
lar para meters (Section~ |3.3| ), and dust-production (Sec- 
tion |3.4|) of the a-AGB stars to the rest of the Magellanic 
Cloud TP-AGB stars. All evidence suggests that the a- 
AGB stars are low-mass stars at the very end of their 
evolution. In Section [4j we compare the results to exist¬ 
ing empirical mass-loss models and to stellar evolution 
models, and we suggest a refinement of the canonical 
TP-AGB classification scheme. 


2. DATA & SOURCE CLASSIFICATION 


Fig. 1.— Color-magnitude diagrams of LMC TP-AGB stars, as 
sele cted in Bll. Uvver Pan el: J — K s versus 7G diagram showing 
the |Cioni et al.| (|2006| and |Blum et al.| (|2006j) divisions between 
C-rich, U-nclL, and ‘extreme' (x-)AUh> stars (black lines). Lower 
Panel: J — [8.0] versus [8.0] diagram. The red lines highlight the 
location of O-rich, C-rich, and a-AGB (originally named aO-AGB 
in Bll) - the a-AGB stars form their own sequence in this diagram 
that is separate from the other AGB stars. In both pan els, the color 
d enote s the density of AGB stars. In the canonical | Cioni et al.| 
(2006) classification scheme, the a-AGB stars overlap primarily 
with O-rich AGB stars, separating only when the 8-fam photometry 
is included. 

In Bll, we put forth three possible explanations for 
the seemingly anomalous J — [8] colors of the aO-AGB 
stars. First, they may be a subset of the O-AGB stars 
that have formed, or are just beginning to form, dust in 


Table [2] summarizes the final source counts used in our 
analysis, classified as described in this section. 

2.1. Photometry 

We use the photometry databases compiled by the Sur¬ 
veying the Agents of Galaxy Evolution (SAGE) Legacy 
programs (Meixner et al. 2006 Gordon et al. 2011), 
which includes photometry from the Magellanic Clouds 
Phot ometric Survey (MCPS, UBV_P Zar itsky et al. 
2002), 2MASS (JiLiC s ;|Skrutskie et al. 2006|), and Rpitzer 
(3.6, 4.5, 5.8, 8.0, and 24 /1 m) photometry. To this, we 
add available archival and literature photometry from 
the Wide-Field I nfrared Sur vey Ex plorer (WISE: 3.4, 4.6, 
11.6, and 22 /im; Wright et al.|2010 ) and the AKARI tar- 
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TABLE 2 

Classification Summary 


a-AGB sample includes O-AGB stars redder in J — [8] 
than the line described by eq. 3 in Bll: 




Photometric Class 



O-AGB 

C-AGB 

x-AGB 

a-AGB 

SMC: 

Total 

2453 

1708 

349 

1266 

O-rich spectrum 




100 

C-rich spectrum 




122 

S-type spectrum 




23 

Unknown spectrum 




28 

No Spectrum 

2453 

1708 

349 

993 

In OGLE-III 

1784 

1459 

229 

1006 

LMC: 





Total 

10 810 

6184 

1389 

6342 

O-rich spectrum 

1041 

44 

4 

419 

C-rich spectrum 

87 

1501 

112 

122 

Unknown spectrum 

20 

23 

142 

7 

No spectrum 

9662 

4616 

1131 

5794 

In OGLE-III 

7445 

4538 

710 

4939 

Note. — These 

are the 

source 

counts and clas- 


sifications included in our analysis. The columns 
represent the photometric classifications, the rows 
indicate spectral classifications and whether a source 
is included in the OGLE-III spatial coverage. The 
LMC numbers exc lude the kinematically-distinct 


population (Section 2.3). 


geted surve ys of the SMC and LMC (3.2, 4.1. 7, 11, 15, 
and 24 /im; Ita et al. 2010 Kato et al. 2012). Table fl] 
summarizes the photometric data used in this study. 

The galaxies were covered by Spitzer twice, 3 months 
apart. The photometry used here is derived by com¬ 
bining these two epochs into a single, co-added stack. 
The magnitudes are therefore a flux average between the 
two epochs, making them less susceptible to stellar pul¬ 
sations. All magnitudes are corrected for foreground ex¬ 
tinction using the same parameters listed in Table 1 from 
Bll and the wave length dependence derived by |Indebe^ 
touw et al.l (|2005h . 


2.1.1. Photometric Classification 

A detailed description of the selection criteria for TP- 
AGB stars is discussed by Bll, and we use identical clas¬ 
sifications in this work. In brief, the NIR J — K$ vs. K$ 
color-magnitude diagram (CMD; Fig. T|) is used to ini¬ 
tially sep arate the C-AGB and O-AGB stars, following 


Cioni et al. (2006). We also classify “extreme” dusty 
AGB stars (x-AGB) as those with J — [8.0] > 3.4 mag or 
J — Ks >2.2 mag. Many studies indicate that the ma¬ 
jority ofthe x-AGB stars in the MC s are C-rich (Woods 


et al. 2011 
Ruffle et a 


I van Loon et al.| 2008| Blum et akT ^OM 
2015]), though obscured O-rich stars can 


also reach very re d colors (e.g., the LMC OH/IR star 
IRAS 05298—6957; |Wood et al.||1992] |van Loon et al. 
2001|). We emphasize that this classification scheme 


includes only stars that are brighter than the 3.6 /am 
TRGB, which is dominated by AGB stars undergoing 
thermal pulsations; we do not include the fainter, less 
evolved early-AGB stars since these are difficult to dis¬ 
tinguish from red giant branch stars. Hereafter, we use 
the term AGB interchangeably with TP-AGB. 

The a-AGB stars are isolated from the initial NIR- 
selected C- and O-AGB populations based on their posi¬ 
tion in the J — [8] vs. [8] CMD (Fig. [lj. Specifically, the 


[8] = A — (11.76 x (J — [8])), 


(i) 


with A = 27.95, and C-AGB stars bluer than the same 
line with A = 31.47. The LMC includes >6000 a-AGB 
stars, and they comprise 26% of the entire LMC TP- 
AGB population. In the SMC, there are ^1200 a-AGB 
stars, and they account for a slightly smaller fraction of 
the SMC’s TP-AGB population (21%). 

These photometric classification are not precise, and 
there is a small amount of cross-contamination between 
each class of stars. We also stress that the a-AGB and 
x-AGB classes are defined purely by the photometric col¬ 
ors; both classes can co ntain both O-rich and C-rich AGB 
stars. In Section 2.3. 1[ we compare the photometric clas¬ 
sifications to spectral classifications. 


2.2. Variability 

To investigate the pulsation properties of the AGB 
stars, we include data from the Optical Gravitational 
Lensing Experiment (OGLE-III; Soszynski et al. 2008) 
surveys of the Magellanic Clouds'! The OGLE survey 
monitored the MCs in the /-band for approximately 8 
years, covered the stellar bars in both galaxies, and de¬ 
rived periods and amplitudes for every variable star. 



Fig. 2.— 3.6 period-luminosity relationship for the LMC, 
showing the source density. Periods are from the OGLE-III survey 
an d magnitudes are from Spitzer. Sequences are labeled follow¬ 
ing [Fraser et al. | (|2005|); we assigned stars to a pulsation sequence 
according to the solid lines. 

The 3.6 /im period-luminosity ( P-L ) relationship for 
the LMC is shown in Figure [2} The sequences first de- 
scribed by | Woo d et al.| (fl999 ) and further subdivid ed by 
Ita et al. (|z(J U4| are clearly evident. We adopt the Fraser 
et al. ~(2O05|) numbering scheme: sequences 1-4 describe 
an evolutionary sequence, with stars evolving to longer 
periods and brighter magnitudes as they age. Stars on 
Sequence 1 are pulsating in the fundamental mode and 
higher overtone pulsators occupy Sequences 2-4. Stars 
on the upper end of Sequence 1 are often called Mira 
variables after the prototype Mira, oCeti. The cause of 
the long secondary periods in stars on Sequence D is un¬ 
known, but might have to do with binar ity, non ra dial p ul- 
sation, or nonspherical geometry (e.g., Wood & Nicholls 














































































4 


Boyer et al. 


TABLE 3 

Spectroscopic Data Summary 


Galaxy 

Instrument 

A 

R 

Date 

SMC 

AAOmega/2dF 

4225 - 4975 A 

3700 

Sep 2011 

SMC 

AAOmega/2dF 

6325-6775A 

8000 

Sep 2011 

LMC 

Hydra-CTIO a 

6850-9150A 

7000 

Nov 2007 


a See Olsen et al. (2011) for a detailed description 
of the Hydra-CTIO observations and data reduction. 


2009). Here, we assign stars to each sequence based on 
the lines drawn in Figure [2j with a shift to fainter abso¬ 
lute 3.6-/im magnitudes in the SMC to account for the 
difference in distance. 


2.3. Spectra 

We obtained medium-resolution optical spectra of 273 
randomly selected a-AGB stars in the SMC barjising the 
AAOmega/2dF multi-object spectrograph (|Lew is et al7 
2002| |Saunders et al.||2004| |Sharp et al.||200b|) on the 
Anglo-Australian Telescope. The total exposure time 
was 3 hr, and reductions were carried out with the 2dF 
pipeline ( 2dFdr , version 4.66, 2011 October). Medium- 
resolution optical spectra for AGB stars in the LMC were 
obtained from 2007 November 20-26 at the CTIO 4 m 
Blanco tel escope with th e Hy dra-CTIO multi-fiber spec¬ 
trograph (Barden & Ingerson 1998). The LMC sample 
comprise s 3791 AGB stars in both the stellar bar and the 
disk. See Olsen et al. (2011) for a deta iled description of 


the observations and data reductions. Olsen et al. (|2011) 
used the LMC spectra to identify a kinematically- and 
chemically-distinct population of stars among the LMC 
sample. They concluded that this distinct population 
was accreted from the SMC during a close passage be¬ 
tween the clouds. We exclude this accreted population 
from the analysis in this paper; 3522 AGB stars, includ¬ 
ing 548 a-AGB stars, remain in our sample (Table [2J. 

Spectroscopic data are summarized in Table [3] fig¬ 
ure |3[ shows the map of the a-AGB spectra in both galax¬ 
ies, classified as described in the next section. 


2.3.1. Spectral Classification 

The AAOmega/2dF spectra (SMC; a-AGB stars only) 
are classified by eye as C (C-rich) or M (O-rich) stars 
based on the presence of C 2 +CN or TiO, respectively 
(Fig. Bl). S-type AGB stars are identified by the pres¬ 
ence 01 ZrO. This classification scheme results in 122 C 
stars, 100 M stars, 23 S stars, and 28 unknown sources 
(Table [2]). We list the SMC S stars in Table [4j 
The Hydra-CTIO spectra (LMC AGB stars only) are 
classified by computing band-strength indices covering 
TiO and CN to measure the strength of these spectral 
features (Fig. [5l Table [5]). We calculate the index for 
each feature asiollows: 


r / 771mm _ rimin \ / 771mm / r \\ 

1 V^cont. ^band//^cont. ’ \^J 

where F con t. and F \> an d are the normalized flux in the 
continuum and band ranges, respectively (Table [5b. The 
only S-type molecular features present in the wavelength 
region covered by the CTIO spectra are due to LaO at 
7385A and 7406A, but these features are weak and can¬ 
not be reliably identified in this sample. We are thus 


LMC 



RA (J2000) 

SMC 



Fig. 3. — Map of the spectrally-classified a-AGB stars in the 
LMC (upper panel) and SM C (lower panel) . The background i s 
the 8-^m Spitzer image (Meixner et al. 2006 Gordon et al. 2011 ). 
The a-AGB stars are marked as blue circles (O-rich), red tri angles 
(C-rich), and yellow squares (S-type), as classified in Sect ion [2. 3. 1| 


limited to classifying the LMC stars as either C-rich 
(C/O > 1) or O-rich (C/O < 1). We classify sources with 
(/cn) > 0.1 as C-rich and sources with (/cn) < 0.05 and 
(Ano) > 0 as O-rich. Figure 6] shows the classification 
of the LMC spectra comp areat o the photometric clas¬ 
sifications from Bll (Sect. |2.1.1 ). The Bll photometric 
classifications agree with the spectral classifications for 
96% of the C-AGB stars and 91% of the O-AGB stars 
(see Table [2|. 

There are 118 spectra that do not meet one of the 
spectral index criteria described above, and we classify 
these by eye (i.e., those in the white areas of Fig. [6] — 
note that some fall outside of the plotted range). We also 
visually inspect and, if necessary, adjust the classification 
for the stars with a spectroscopic index that disagrees 
with the photometric classification (see Fig. [6]), which 
included 55 and 206 photometrically classified^ O-AGB 
and C-AGB stars, respectively. 

The x-AGB stars are faint in the optical, resulting 
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7000 7500 8000 8500 X [A] 9000 


Fig. 4.— Example AAOmega/2dF specta of SMC a-AGB stars 
in normalized flux, binned to a resolution of R ~ 400. Molecular 
features used to classif y stars are marke d. C? and CN bandheads 
are f rom [D avis (198 71, TiO from IV alenti et al.l fl998), and ZrO 
from Garcia-Herna ncfez et al.j ( |200Yi r 


Fig. 5.— Example CTIO specta (in normalized flux) of a C 
a nd M AGB star in the LMC, binned to R ~ 400 (jOlsen et al. 
|2011|. Molecular fe atures used to classify stars are marked. (J IN 
bandheads are from | Davis |( |1987| >, TiO from | Valenti et al~] ( |1998| >. 


TABLE 4 

Spectroscopically-identified S-type Stars in the SMC 



Designation 

(SSTISAGEMA) 

[3.6] 

(mag) 

[8.0] 

(mag) 

Period 

(d) 

1 

J004234.17-725300.6 

11.294=0.01 

10.92i0.03 


2 

J004428.44-730927.9 

11.78T0.02 

11.70i0.02 

80 

3 

J004526.04-721432.9 

ll.66i0.03 

11.39i0.03 


4 

J004547.21-730535.0 

11.31i0.01 

11.18i0.01 

110 

5 

J004618.11-732652.4 

12.44i0.02 

12.20i0.05 

401 

6 

J004712.99-724637.6 

11.38i0.01 

11.03i0.02 

653 

7 

J004721.34-724834.3 

11.65i0.01 

11.45i0.02 

91 

8 

J004731.02-732942.6 

11 .20i0.01 

11.04i0.02 

165 

9 

J004745.83-732719.6 

11.51i0.01 

11.33i0.02 

95 

10 

J004755.46-732855.3 

11.42i0.01 

11.36i0.02 

82 

11 

J004808.00-724406.3 

11.65i0.01 

11.45=b0.02 

685 

12 

J004833.34-725953.1 

11.62i0.02 

11.38=b0.02 

88 

13 

J004955.64-722755.3 

11.70i0.03 

11.49=b0.04 

89 

14 

J005005.38-730500.1 

11.85i0.01 

11.65=b0.02 

116 

15 

J005342.03-720007.5 

11.86i0.03 

11.73±0.04 

99 

16 

J005404.72-723727.8 

11.37i0.01 

11.27±0.02 

89 

17 

J005418.85-725905.5 

11.37i0.02 

11.25±0.02 

105 

18 

J005458.13-725233.2 

11.90i0.02 

11.67±0.02 

659 

19 

J005510.79-721734.0 

11.44i0.01 

11.33±0.02 

109 

20 

J005744.57-721509.1 

10 .66i0.01 

10.52=b0.02 

123 

21 

J005746.46-723130.4 

11.76i0.01 

11.58±0.02 

93 

22 

J010119.74-725149.3 

11.62i0.02 

11.44±0.02 

606 

23 

J010133.86-732135.8 

11.61i0.03 

11.41±0.04 

92 


Note. — The full photometric catalog is available on¬ 
line, including UBVI from MCPS, JHK S from 2MASR 
and IRAC and MIPS magnitudes from SAGE (see Table [l] 
for references). The OGLE-III periods and amplitudes are 
also included in the online version. 


in noisy spectra that make spectral index analysis dif¬ 
ficult. We therefore check by eye, and revise if nec¬ 
essary, the spectral classification of all x-AGB stars. 
We are able to identify carbon or oxygen molecular 
features in 116 of the 258 x-AGB stars in the sam¬ 
ple. All but four of these 116 are C-rich; the O-rich x- 
AGB stars include SSTISAGEMA J045414.27-684413.9, 
SSTISAGEMA J051333.64-683644.4, SSTISAGEMA 
J052710.19—693626.9 (WOHG339), and SSTISAGEMA 


TABLE 5 

Spectral Indices for LMC Spectra 


Feature 

Band 

A (A) 

Continuum 

A (A) 

TiO-1 

7051—7057 

7014—7051 

TiO-2 

7122—7128 

7014—7051 

TiO-3 

7700—7750 

7480—7530 

CN-1 

7250—7300 

7215—7240 

CN-2 

7860—7920 

7800—7845 


J044525.22—704648.3. The first of these is a known O- 
rich Mira variable, with strong 10-/im silicate emission 
(alternate identifi cations: LI-LMC153, SP7730-6, and 
IRAS 04544-6849; Loup et al.|1997||van Loon et al.|1999 

".. TG3~" ’ ’ 


Dijkstra et al. 12005 ). WOHG 339 is also known to show 

weak silicate emission ( Sloan et al.||2008 ). 


3. ANALYSIS: PROPERTIES OF A-AGB STARS 

In this section, we summarize the observable properties 
of a-AGB stars. 


3.1. Chemistry 

Table [2] summarizes the spectral classifications of the 
a-AGB stars. In the LMC, 76% of the a-AGB spectra are 
O-rich, and this fraction decreases to 45% in the SMC 
(including M and S stars). In Figure [7| we show the 
histogram of a-AGB stars in J — [8]. In both galaxies, 
a-AGB stars are equally likely to be C-rich or O-rich at 
red colors (J — [8] > 1.75 mag in the LMC and >1.55 mag 
in the SMC), and O-rich stars dominate at bluer colors, 
though to a lesser degree in the SMC. 

The high fraction of C-rich stars among the a-AGB 
sample indicates that caution is necessary when using the 
standard J — K$ color method to separate O-rich and C- 
rich stars. Although the spectra indicate that only about 
half of the a-AGB stars are O-rich in the SMC, more than 
96% of them would classify as O -rich with the color cuts 
suggested by Cioni et al. ( 2006). The exact placement of 
these cuts also varies depending on the chosen foreground 
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Fig. 6. — TIO vs. CN spectral indices for LMC stars compared to 
the photometric classifications from Bll. Each subpanel includes 
sources that were photometrically classified as indicated by the 
label in the upper right. The blue and red shaded regions illustrate 
the spectral index classification employed here for O-rich and C- 
rich sources, respectively. We classified stars in the white area by 
eye. This diagram demonstrates the cross-contamination between 
the photometric C-AGB and O-AGB groups (at the ?^4%-9% level), 
and shows that the a-AGB photometric class includes a sizeable 
fraction of both O-rich and C-rich sources. Also see Table l2l 



Fig. 7.— J — [8] histogram of a-AGB stars in the LMC and 
SMC. Stars classified as O-rich or C-rich by their optical spectra 
are plotted in blue and red, respectively. In the SMC, the blue line 
includes S-type and M-type stars, and S-type stars are also plotted 
separately in orange. S-type stars cannot be identified in the LMC 
sample due to limited wavelength coverage. In both galaxies, the 
bluest a-AGB stars tend to be O-rich. 


extinction correction and the assumed metallicity of the 
galaxy, and slight shifts result in potentially large uncer¬ 
tainties in the C/M ratio. By assuming that all a-AGB 
stars are O-rich, Bll found C/M ratios of 0.56 in the 
SMC and 0.42 in the LMC. After correcting for the frac¬ 
tion of a-AGB stars that are spectrally classified here as 
C-rich, we find revised C/M ratios of 0.87 and 0.56, re¬ 
spectively. Relying on the J — K$ color is therefore not 
ideal for separating C and M stars. 

Between M and C stars lie the S-type stars; as an AGB 
star evolves, carbon and 5-process elements are dredged 
up to the surface. As C/O in the stellar envelope ap¬ 


proaches unity, TiO bands in M stars begin to disappear 
in favor of ZrO bands, which are prominent in S-type 
AGB stars (though C/O can be near unity even in the 
absence of ZrO bands). As more carbon is dredged up, 
free oxygen is tied up into CO and carbon molecules dom¬ 
inate the stellar spectrum. ZrO features are prominent 
in two types of S stars: intrinsic and extrinsic. Intrin¬ 
sic S stars follow the evolution sequence just mentioned. 
Extrinsic S stars are enriched in carbon and 5-process el¬ 
ements through accretion from a binary companion that 
has already finished its AGB evolution. Extrinsic S stars 
are less evolved than intrinsic S stars, and are usually 
found on the RGB or among early-AGB stars that have 
not yet begun to thermally pulse. To identify intrinsic 
S stars, studies rely on (1) the presence of technetium, 
which will have decayed in the time since the mass trans¬ 


fer epis ode that created the extrinsic S star (e.g., Joris- 


2003), (2) infrared colors that indicate the presence of 


dust, whic h cannot form in less-evolved RGB and E-AGB 
stars (e.g., |Groenewegen 1993 Jorissen et al.||1993 Otto 


et al. 


2011 ) 


Groenewegen 1993| (Jorissen et al.||1993 

), and (3) high luminosities that are c 


aarac- 


teristie^Tthermally-pulsing AGB stars (e.g., Guandalini 
fc Busso||2008 ). 


Yang et al. (|2006) identified Galactic intrinsic and ex- 
trinsic S-type AGB stars based on their IR colors. Bll 
noted that the_colors ofjdie a-AGB stars are consistent 
with the Yang et al. (2006) extrinsic S stars, but with¬ 
out spectroscopy, we could not confirm or rule out the 
presence of S stars among the a-AGB sample. Here, we 
find ZrO features in the spectra of only 9% of our SMC 
sample of a-AGB stars, so an S-type nature cannot ex¬ 
plain the red J — [8] colors of the entire a-AGB popula¬ 
tion. Moreover, despite their similar colors to Galactic 
extrinsic S stars, the S stars detected here a re lik ely to be 
intrinsic because their lum inos ities (Section pL3| and pul¬ 
sation properties (Section |3.2| ) are consistent with AGB 
stars undergoing thermal pulses. Th e overl ap in color be¬ 
tween the a-AGB stars and the Yang et al. (2006) extrin¬ 
sic sample suggests that IR colors may not be a reliable 
way to separate intrinsic and extrinsic S stars. 


3.2. Pulsation Properties 

Figure M shows the LMC 3.6 /am P-L diagram for each 
type of (photometrically-classified) AGB star; the pro¬ 
portions of stars on each sequence in the SMC are sim¬ 
ilar. Table [6] shows that the a-AGB stars preferentially 
occupy Sequences 1, 2, and D (we discuss Sequence D 
stars in Section (4~4| ) . Most a-AGB stars are classified by 
OGLE as semi-regular variables (SRVs), which have pul¬ 
sation properties similar to Miras, but with less regular 
light curves and smaller E-band amplitudes. 

In both galaxies, the majority of the fundamental¬ 
mode pulsators (Sequence 1; those at the end of the 
TP-AGB phase) are photometrically classified as C-rich 
(either C-AGB or x-AGB). Only 9-10% of them are pho¬ 
tometrically classified as O-rich, and 15-30% are a-AGB 
stars. Almost 90% of the O-AGB stars have shorter pe¬ 
riods (Sequences 3 and 4), indicating that the O-AGB 
class as defined in Bll comprises mostly stars at the be¬ 
ginning of the TP-AGB phase. 

In Figure [9| we show the 3.6 /am P-L relationship for 
a-AGB stars classified via their optical spectra. In both 
galaxies, the O-rich a-AGB stars (with TiO or ZrO fea¬ 
tures) lie on sequences 2 and D. In the LMC, they also 
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Fig. 8 .— 3.6 /im period-luminosity diagrams for the LMC (greyscale), separated by photometrically-classified spectral type (black dots). 
For the O-AGB and C-AGB stars, we exclude sources that are spectrally classified as the other. The a-AGB stars are shown in Figure [9] 


TABLE 6 

3.6 fiM P-L Relationship Statistics 


AGB 

class 

Seq. 4 

N (%) 

Seq. 3 

N (%) 

Seq. 2 

N (%) 

Seq. 1 

N (%) 

Seq. D 

N (%) 

a 

213(4) 

Photometric Classifications - 
LMC 

448(9) 2135(43) 1140(23) 

1003(20) 

O 

2756 (37) 

2248 (30) 

1066(14) 

373 (5) 

1002(13) 

C 

41(1) 

297 (7) 

2096(46) 

1634(36) 

470(10) 

X 

9(1) 

10 (1) 

23(3) 

656(92) 

12 (2) 

a 

56 (6) 

118(12) 

SMC 

317(32) 

160(16) 

355(35) 

O 

566(32) 

443(25) 

224(13) 

103 (6) 

448(25) 

C 

7(1) 

53(4) 

588 (40) 

603(41) 

208(14) 

X 

0 (0) 

2 (1) 

4(2) 

219(96) 

2 (1) 


- Spectral Classifications; a-AGB only - 
LMC 


o 

6(3) 

15(6) 

99 (43) 

63 (27) 

48(21) 

c 

7(9) 

19(25) 

19 (25) 
SMC 

11(15) 

19(25) 

o 

5(6) 

5(6) 

29 (32) 

9(10) 

42 (47) 

c 

14(13) 

28 (26) 

24 (22) 

13(12) 

29 (27) 

s 

0 (0) 

1 (6) 

12 (57) 

3(14) 

5(24) 


Note. — Position of AGB stars on the 3.6 /am P- 
L diagram. The percentage of the total population of a 
particular subtype is shown in parentheses. For example, 
213 of 4939 a-AGB stars (4%) are found on Sequence 4 in 
the LMC (see Table [2|. For an independent acc ounting 
of the pulsation sequences, see|Spano et al. (2011). 


occupy Sequence 1. Pulsation models (Wood 2015) show 
that stars move to longer periods as their luminosities 
increase, suggesting that many high-L a-AGB stars on 
Sequence 2 will evolve into C stars as they cross to Se¬ 
quence 1 and that most of the O-rich a-AGB stars on 
Sequence 1 will remain O-rich. The absence of a-AGB 
stars on Sequence 1 in the SMC reflects the efficiency of 
this transition in metal-poor environments. 

We have identified 24 S-type stars in the SMC, and 
22 of them are classified as AGB stars by the OGLE-III 
survey (the remaining 2 fall just outside the OGLE-III 
spatial coverage). The OGLE light curves thus provide 
additional evidence that the 5-process elements evident 
in their spectra are from the dredge-up of carbon during 
a thermal pulse (intrinsic S-type), and not due to mass 
transfer (extrinsic S-type). The clustering of these intrin¬ 
sic S stars on the 3.6-/im P-L diagram is an intriguing 
clue about the dredge-up process. However, the SMC 
spectral sample does not include stars from the C-AGB 
and O-AGB photometric categories, so the frequency of 
S-type stars among those groups remains unknown, as 



Fig. 9.— 3.6 /im period-luminosity diagrams for a-AGB stars. 
Grey dots show all a-AGB stars, and colored symbols mark those 
with spectral classifications. Note the change in y-axis from Fig¬ 
ure [8] and the y-axis change between the LMC and SMC. 

does the position of any such stars on the P-L diagram. 


3.3. Fundamental Parameters: L, T e q, M, and [Fe/H] 

To estimate the luminosities (L) and effective tem¬ 
peratures (T e ff), we fit the spectral energy distribu¬ 
tions (SEDs) of the AGB stars using the photometry 
list ed in Table [l] and following the procedure outlined 
by |McDonald et al. (2009 2012). Each SED was com¬ 
pared to a grid of b t-settl stellar atmosphere models 
(Allard et al. 2011) using y 2 minimization and scaled 
in flux to derive a bolometric luminosity. We assume 
[Fe/H] = —0.3(—0.5), [a/Fe] = +0.3(+0.3), d = 51 kpc 
(61 kpc), and E(B — V ) = 0.085 mag (0.070 mag) for the 
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Fig. 10.— Hertzprung-Russell diagram for stars on pulsation 
Sequences 1 and 2. The x-AGB stars are excluded because circum- 
stellar extinction results in unreliable SED fits. The a-AGB stars 
are intermediate to the O-AGB and C-AGB stars in both T e ff and 
L. 

LMC (SMC). 

Uncertainties are large for stars with significant cir- 
cumstellar dust since this method compares SEDs to 
dustless photospheres. In practice, we find that SED 
fits with temperatures less than 2800 K, encompassing 
most of the x-AGB stars, are unreliable. We therefore 
exclude most x-AGB stars when discussing effective tem¬ 
peratures, luminosities, and pulsation masses. 

We show the Hertzprung-Russell diagram in Figure [lOj 
To compare stars in a similar evolutionary phase, we in¬ 
clude only stars on pulsation Sequences 1 and 2. In both 
galaxies, the a-AGB stars have temperatures and lumi¬ 
nosities intermediate to the O-AGB and C-AGB stars. 
This suggests that a-AGB stars are those that are near 
a transition from O- to C-rich. The temperatures and 
luminosities of the a-AGB stars here are similar to those 
determined for the C-rich and S-type stars in the Hijy 


yarco s catalog using the same method (McDonald et al. 
2012 ). 


To derive the stellar masses, we rely on th e period - 
mass-radius relationship described by Vassiliadis fc| 


Wood (|1993|). This relationship is valid only for stars 
that lie on pulsation Sequence 1 and are therefore pul¬ 
sating in the fundamental mode: 

logP(d) = -2.07+1.94log R/R q —0.9logM p /M 0 , (3) 

with R derived using L rsj The mass derived 

from this relationship is the pulsation mass (M p ), which 
reflects the current mass of the star, so small masses 
could either indicate a low-mass star that has undergone 



1.5 2.0 

Pulsation Mass (M 0 ) 

Fig. 11.— Distribution in pulsation mass for stars on Sequence 1 
in the LMC. Stars included here are photometrically classified. The 
thick gray line shows the sum of all three components, the dotted 
line marks the mean, and the dashed line marks the median. 

little mass loss or a high-mass star that has undergone 
substantial mass loss. 

Figure [TT] shows the distribution in pulsation mass for 
the fundamental-mode pulsators, divided by stellar type. 
Most stars (95%) h ave pu l satio n masses of 0.5-2.5 M®, 
in agreement with Wood (2015), and the a-AGB stars 
show the lowest median pulsation masses. Since the a- 
AG B ma ss-loss rates are not particularly high (see Sec¬ 
tion [+4 ), we can infer that the a-AGB stars have the 


lowest initial masses among the stars pulsating in the 

fundamental mode. _ 

Meta llicitie s were derived using the Ca triplet by I Olsen 
et al.| (|2011), including 36 a-AGB stars and 268 O- 
AGB stars m the LMC. Abundances were not reliably 
measured for any C-AGB or x-AGB stars. The me¬ 
dian metallicity of the O-AGB population is [Fe/H] = 
—0.82 dex, with dispersion (median absolute deviation) 
of 0.38 dex and a median uncertainty of 0.17 dex. The 
a-AGB population is more metal-poor, with a median 
[Fe/H] = —1.05 dex, dispersion 0.36 dex, and median un¬ 
certainty 0.19 dex. The lower metallicities suggest that 
the a-AGB population is older, having formed earlier in 
the LMC’s chemical evolution history. This finding is 
consistent with the lower masses inferred from the pul¬ 
sation. Alternatively, the a-AGB stars may have been 
accreted from the SMC or formed from gas accreted from 
the SMC. 

3.4. Mass Loss and Dust Production 

Fourteen of the LMC a-AGB stars were observed with 
Spitzer’s infrared spectrograph (IRS) from 5-16 /am as 
part of the SAGE-Spec program (Kemper et al. 2010). 
The equivalent pr ogram for the S1V1G does not include 
any a-AGB stars (Ruffle e t al. 20151) . Of the 14 LMC 
stars, four are classified by | Woods et al.| (|2011|) as Stars 
(their AGB nature is unconfirmed) ana the remaining 
ten as O-rich AGB stars, nine of wh ich s how strong sil¬ 
icate emission at 10 /am (see Section 4.6). With only 14 


examples, we cannot draw any broad conclusions about 
the IR spectral features in a-AGB stars, but the presence 
of a silicate feature in these examples demonstrates that 
a-AGB stars are capable of significant dust production. 
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J-[8.0] (mag) 

Fig. 12 .— The color-DPR relationship for LMC stars on Se¬ 
quence 1. The blue and orange points were classified by their spec¬ 
tra, and the black dots are the photometrically-identified a-AGB 
stars (all Sequences). In the LMC, most of the a-AGB stars follow 
the low-mass end of the O-rich branch. 


Riebel et al. (2012) computed the mass-loss and dust- 
production rates (DPRs) for the LMC AGB stars by 
fitting the full optical to mid-IR SEDs to the Grid of 
AGB and RSG Model S (GRAMS; |Sargent et al. 2011 


Srinivasan et al 


l20TTl). 

liar pho 


The GRAMS grid inclu des the 


COMARCS stellar photospheres for C stars (Gautschy- 
Loidl et al.|2004|) and PHOENIX models for O-rich stars 


(|Hauschildt et al. 1999). Silicate dust or a mixture of 
amorphous carbon-f SiC were added to create a grid that 
reproduced the full luminosity and infrared color distri¬ 
bution of the LMC AGB population. We note that lumi¬ 
nosities derived with GRAMS using a simple trapezoidal 
integration are systematic ally lower than those derived 
by the SED fitting (Section [3.3| ) by about 16% for O-rich 
stars and 18% for C-rich stars. 

Fo r a s ubse t of so urces, the GRAMS classification from 


Riebel et al. (2012) disagrees with our spectral classifica¬ 
tion. This includes 29 of the spectrally-confirmed O-rich 
sources that were classified by GRAMS as C-rich, and 
241 of the spectrally-confirmed C-rich sources that were 
likewise classified as O-rich by GRAMS. GRAMS pro¬ 
vides a best-fit O-rich and C-rich model for each star, 
and the GRAMS classification is assigned depending on 
which provides a smaller y 2 , so many misclassifications 
are the result of only a very slight preference for one fit 
over the other. Here, we assign the best-fit C- or O-rich 
GRAMS model that matches the spectral classification. 

The GRAMS grid assumes a wind expansion velocity 
of lAxp = 10 kms -1 , a spherical dust shell geometry, and 
a constant DPR throughout the AGB lifetime. The DPR 
scales linearly with t> exp , which itself scales with luminos¬ 
ity and the gas-to-dust ratio (0) as f exp oc L 0,25 ^ -0,5 . To 
scale the GRAMS DPRs, we use different gas-to-dust ra¬ 
tios depending on the stellar chemistry. For O-rich stars, 
we assume that ip scales with metallicity, so ^lmc h = 400 
for Zlmc = 0.5 Zq and -0© = 200 (here we use the typ¬ 
ically assumed LMC metallicity). For a-AGB stars, we 
scale ip according to the difference in metallicities be¬ 


tween O-AGB and a-AGB stars discussed in Section 3.3 


so V^lmc B = 080. Since C stars make their own carbon, 


TMC 

we expect their gas-to-dust ratios to be near solar, so we 
Tmc 1 = 200. 


use 


Estimated DPRs from SED fitting are highly sensitive 
to the choice of dust optical constants and dust co mpo¬ 
sitions, which may differ at these me tallicities (e.g., Mc¬ 
Donald et al.|2010 ; Jones et al.|2014 ). We therefore note 
that while the relative GRAlVlS-derived DPRs for differ¬ 
ent stellar types ar e valid, absolute D P Rs are inaccurate 
by a factor of ^4 ([Gr oenewege n et al. 2009; Srinivasan 
et al.||2011l Ibargent et al.||2Ull| |lVlcDonald et al.||201l|). 

The DPRs indicate that the a-AGB st ars "are mod¬ 

erate ly strong dust producers (also see [Boyer et al. 
2012). The median DPR of the a-AGB population is 
4 x 10 -11 M e yr -1 (dispersion 2 x 10 -11 M© yr _1 ). This 
is an order of magnitude higher than the median DPR of 
the O-AGB population, which includes many non-dusty 
stars, and a few bright very dusty stars. The C-AGB and 
x-AGB populations are very dusty, with median DPRs 
that are ~2 x and 2 orders of magnitude higher than the 
a-AGB median DPR, respectively. 

Figure 12 shows how the J — [8.0] color varies with 
the DPR; colors represent the spectral classification and 
symbol sizes represent the pulsation mass (only Se¬ 
quence 1 stars are plotted, so the low-DPR O-AGB pop¬ 
ulation, which resides on Sequences 3 and 4 are not in¬ 
cluded). The relationship between J — [8] and DPR is 
bimodal, with O-rich stars reaching higher DPRs at bluer 
colors than C stars and a clean separation between the 
two spectral types. O-rich stars increase in mass as the 
DPR increases, but the masses of C-rich stars is mixed 
along the entire range of DPRs. Fits to each branch 
shown in Figure [12] are of the form: 


logD = 


x 

(J — [8.0]) + y 


(4) 


for stars with J — [8.0] > 1 mag. For C-rich stars, x = 
—3.60, y = —0.58, and z — —8.39. For O-rich stars, 
x = —3.54, y = —0.42, and 2 = —7.30, though it is 
unclear whether this relationship applies to O-rich stars 
with J — [8.0] > 3 mag, since we have few examples of 
O-rich stars at those extreme colors. 

Overplotted as dots in Figure [T2] is the entire popula¬ 
tion of a-AGB stars, including all pulsation sequences. 
The bulk of the a-AGB stars smoothly continue the O- 
rich branch to low DPRs and masses, indicating that the 
a-AGB stars are a subset of the general (dust-producing) 
O-rich AGB population in the LMC. 


4. DISCUSSION 

The properties of the a-AGB stars presented in the 
previous section all suggest that they are low-mass, very 
evolved dusty AGB stars. The broad spatial distribu¬ 
tions of the a-AGB stars compared to other AGB stars 
in the SMC presented in Bll also suggests that a-AGB 
stars originate from an older stellar population. The star- 
formation histories of the Magellanic Clouds have been 
computed globally by Harris fe Zarit s ky (120041 2009), 
and in smaller fie lds byjRiibele et al. (|2U12 |2UloJ) and 
Weisz et al. (2013). While there is strong variation m the 
star-formation histories regionally, global histories show 
a peak near £ age = 2-3 Gyr, which corresponds to an ini¬ 
tial stellar mass of 1-1.3 M©. The a-AGB stars likely 
originated from this star-formation event. This mass 
range encompasses the lower-mass limit for st ars t hat 
will become C-rich during their evolution (e.g., Marigo 
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Fig. 13.— The relationship between pulsation period and mass 
loss for AGB stars pulsating in the fundamental mode (Sequence 
1 ) in the LMC. The upper panel shows C-rich stars and the lower 
panel shows O-rich stars. Large symbols/darker colors indicate 
larger pulsation masses for C-AGB, O-AGB, and a-AGB stars; x- 
AGB stars are plotted with the same colors/sizes since their masses 
could not be reliable estimated. In both panels, we include stars 
classified via photometry, and use the optical spectra to eliminate 
interlopers. The straight lines are from|Vassiliadis & Wood] (1993)) 
for stars with pulsation masses < 2.5 Mq (dashed line) and = 
5 Mq (dot-dashed line), derived by fitting Galact ic C-rich and O - 
rich Miras. The curved thin line is from|Groenewegen et al.|(fl998), 
derived from Galactic C-rich Miras and scaled down to account tor 
the difference in optical constants (see |Srinivasan et al.|20 11). 


et al. 2013). 


in this section, we explore additional evidence support¬ 
ing this low-mass hypothesis, compare the characteristics 
of the a-AGB stars to their more massive counterparts, 
and discuss the opportunities that arise with a simple 
photometric technique for distinguishing low-mass TP- 
AGB stars. We note that the O-AGB population, as de¬ 
fined in Section [2. 1.1 1 can be separated into more evolved 
(pulsation Sequences 1 and 2), and less-evolved (pulsa¬ 
tion Sequences 3 and 4) TP-AGB stars. The following 
sections exclude the less evolved O-AGB stars, except 
where explicitly stated. 


4.1. Dust Production and Pulsation 

Several works point to pulsation as one of the drivers of 
AGB dust production and mass loss, and we now explore 
the influence of pulsation on the dust-production rates of 
a-AGB stars. The influence of pulsation is clear: the me¬ 


dian DPR increases by an order of magnitude from Se¬ 
quences 4 to 3, and again from Sequences 3 to 2, only sur¬ 
passing 10“ 11 M Q yr _1 on Sequences 2, 1, and D, where 
the a-AGB stars reside. 

Vassiliadis & Wood ( 1993[ VW93) found that for 


Galactic G-rich and O-rich AGB stars pulsating in the 
fundamental mode, the mass-loss rate increases exponen¬ 
tially with period as: 


log M{Mq yr -1 ) = -11.4 + 0.0123P(days) (5) 

until P = 500 days, when the star enters the superwind 
phase and the mass-loss rate reaches a maximum near 
10 -4,2 M 0 yr -1 . For stars more massive than 2.5 M©, 
their models suggest that this relation changes to 


log M = -11.4 + 0.0125[P - 100(M/M q - 2.5)]. (6) 

The VW93 parameterization is shown for the LMC in 
Figure M alon g with the scaled dust-production rates 
from Riebel et al.| (2012) and the periods from OGLE, 
using the spectral classifications to exclude stars that 
are photometrically misclassified. We convert the gas 
mass-loss rates from VW93 to DPRs using a solar gas- 
to-dust ratio (^ = 200) for C stars and scaling the ratio 
with metallicity for O-rich stars (^ = 400). Equation [5] 
provides a reasonably good description of the behavior of 
C-rich stars. The predicted trend for more massive stars 
to move to longer periods is evident in Figure p~3] for the 
O-rich stars. _ 

Groenewegen et al. (1998) (and later Groenewegen 


et al. It r fit the spectral energy distributions of several 


Galactic C-rich Miras to revise the VW93 M-P relation¬ 
ship. They find 


log M = (4.08 ±0.41) log P - (16.54 d= 1.10), (7) 


which we show in Figure [l3j scaled down by a factor 
of 4 to acco unt for differences in op tical constants as 
described by Srinivasan et al. (2011) and by the gas- 
to-dust ratio {ifo = 200/ ^o = 400). The LMC AGB 
stars do not deviate strongly from the G98 Galactic M- 
P relationship, despite their lower metallicities (Zlmc ^ 
0.5Z o ). 


The LMC a-AGB stars on pulsation Sequence 1 (and 
thus included in Fig. [l3|) are primarily O-rich. Their po¬ 
sition to the left of theaashed line from VW93 indicates 
masses <2.5 M 0 (the O-AGB stars within the a-AGB 
locus could easily also be classified as a-AGB stars with 
a slight shift in the color classifications - see Section |4.5| ). 
Unlike the more massive O-rich stars and the x-AGS/G- 
AGB stars, the a-AGB DPRs are not strongly dependent 
on the pulsation period: their DPRs span two orders of 
magnitude over the full a-AGB period distribution. The 
DPRs of the x-AGB and massive O-AGB stars, on the 
other hand, show a strong positive correlation with both 
pulsation periods and amplitudes. These a-AGB star 
characteristics suggest that once a low-mass star is pul¬ 
sating in the fundamental mode, the pulsation plays only 
a small role in its dust production, at least compared to 
their more massive counterparts. 


4.2. Comparison to Stellar Evolution Models 
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Fig. 14.— 24 fi m CMD for AGB stars on Sequences 1 and 
2. Larger/darker symbols indicate a higher bolometric luminos¬ 
ity. Plotted stars are classified photometrically, using the available 
spectra to mini mize m isclassifications. Evolutionary tracks from 
|DeirAgli et al.| 12015]) at Z = 0.004 (SMC) and 0.008 (LMC) are 
overplotted tor the indicated initial stellar masses. The LMC a- 
AGB stars correspond to <1.5 Mq tracks, while in the SMC, they 
must have initial masses <1.25 Mq, and/or rapidly transition to 
a C-rich atmosphere. The region to the ri ght of t h e da shed line is 
occupied by the “F” finger, noted by |DeirAgli et al. (2015]) and 
encompasses the dustiest a-AGB stars. 


Dell’Agli et al. (2015|) (hereafter D15) computed theo- 
retical evolutionary tracks of AGB stars, including dust 
formation. They found that the masses of AGB stars can 
be inferred by their positions in a Syitzer [8.01 — [24] vs. 
[24] CMD. In the upper panel of Figure 14, we show the 
CMD of the LMC for Sequence 1 and zs tars. As stars 
increase in bolometric luminosity (Section |3.3| ), the size 
of the symbol increases and the color darkens. The most 
luminous (and by extension, the most massive) stars are 
found at the bright/red end of the CMD, following the 
predictions of D15. 

The a-AGB stars occupy a finger at faint 24-/im mag¬ 
nitudes, with a broad range in [8.0] — [24] color. In D15, 
this region of the CMD is populated by stars with initial 
masses <1.5 M 0 , prior to transitioning from O-rich to 
C-rich, which they designate “F” stars (dashed line in 
Fig. fl4k). Their models indicate that “F” stars evolve 
along this finger to brighter 24 /im magnitudes before 
shifting horizontally to the C star sequence when dredge- 
up causes C/O to exceed unity at the very end of the 
AGB phase. Stars with lower masses, and thus insuf¬ 
ficient dredge up for C/O > 1, would remain on this 
finger. Indeed, dusty AGB stars in Globular Clusters 
( M\ < 1 Mq) also occupy the “F” region (Boyer et al. 
2008 2009 McDonald et al. 2012). The D15 models 
support our assertion that the a-AGB are low mass stars 
with insufficient dredge-up: those with the lowest masses 
will remain O-rich, while those with higher masses will 
produce significant quantities of O-rich dust prior to be¬ 
coming C stars. 

At lower metallicities, the lower-mass limit for the 
O- to C-rich transition decreases, and the D15 models 
predict that the F sequence should therefore be more 
sparsely populated in metal-poor environments. Indeed, 
the “F” finger is not pronounced in the more metal-poor 
SMC (Fig. 14). This is partially caused by limited sen¬ 
sitivity at 2M /im, which is 0.4 mag worse in the more 
distant SMC. However, since the a-AGB stars increase 
in brightness at 24 /im as color increases, any in the “F” 
region of the CMD in the SMC are detectable. We thus 
conclude that metallicity is the dominant cause of the 
empty “F” region in the SMC. This effect may be com¬ 
pounded by less-efficient dust production in O-rich stars 
at low metallicity due to a lack of condensation seeds, 
which will result in fewer O-rich stars with 24-/im ex¬ 
cess. Based on the Z = 0.004 D15 models, any low-mass 
stars in the SMC that will remain O-rich must have initial 
masses < 1.25 M 0 , in good agreement with the masses 
inferred from the star-formation histo ries. 

The Padova TP-AGB model^] (Marigo et al. 2008) 
generally follow the D15 models, with slightly higher 
mass transitions. Here, we use the option with 60% sil¬ 
icate + 40% A10 x for O-rich stars and 85% amo rphous 
carbon + 15% SiC for C-rich stars (Groenewegen 2006). 
These models show that, in the LMC, all stars with ini¬ 
tial mass <1.5 M 0 remain on the “F” branch marked in 
Figure 14 transitions over to the C-star side (left side) 
of the diagram occur only for stars with >1.6M 0 . The 
Padova models indicate that this transition happens at 
Mi « 1.3 Mq in the SMC. 

4.3. Probing the lower mass bound for C stars 


1 http://stev.oapd. inaf.it/cgi-bin/cmd 
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Fig. 15. — Luminosity functions for a-AGB stars and 
photometrically-classified C-rich stars (including C-AGB and x- 
AGB stars). We used the spectral classifications to minimize con¬ 
tamination from O-rich sources in the C-rich luminosity function. 
The increase in C-rich stars and decrease in a-AGB stars from Se¬ 
quence 2 to 1 over the same luminosity range suggests that some 
Sequence 1 C-AGB stars originate from the Sequence 2 a-AGB 
sample. 


The a-AGB photometric classification appears to in¬ 
clude initial stellar masses that straddle the lower-mass 
limit for the O-to-C transition caused by the dredge-up 
of carbon. They therefore include both stars that will 
undergo this transition and stars that will remain O- 
rich. We can probe this transition several ways using the 
a-AGB stars. 

First, this bounda ry is reflected in the pulsation masses 
estimated in Section [3731 for stars pulsating in the funda¬ 
mental mode (Sequence 1). The median pulsation mass 
for a-AGB stars is 1.14 M 0 in the LMC and 0.94 M 0 in 
the SMC, with disp ersions of 0.21 and 0.18 M 0 , respec¬ 
tively (Section 3.3). 

Second, the luminosities of S-type stars can probe the 
O-to-C mass transition. In Figure [9| we show the spec¬ 
tral classifications of a-AGB stars on the 3.6-/im P-L 
diagram. In the SMC, we see a curious pile-up of S- 
type a-AGB stars on Sequence 2 near [3.6] = 11.3- 
11.8 mag, corresponding to log(L/L 0 ) = 3.6-3.7 and 
—7.6 < M 3i6 < —7.1 mag. Ass uming these are all in¬ 
trinsic S-type stars (Section |3.2| ), it suggests that stars 
fainter than this cutoff on Sequence 2 will likely remain 
O-rich as they evolve to Sequence 1. Indeed, this is sup¬ 
ported by the relative lack of faint C-rich a-AGB stars 
on Sequence 1 in both galaxies. This scenario would be 
confirmed if LMC S-type stars showed a similar S-type 
star limit, but at a brighter luminosity. Unfortunately, 
the only s-process molecular feature covered by the LMC 
spectra is LaO, which is too weak for reliable S-type clas¬ 
sification. 

Third, a comparison of the luminosity functions of a- 
AGB stars and C-rich AGB stars on different pulsation 
sequences could be used to calibrate the O-to-C transi¬ 
tion in stellar evolution models. We show these functions 
in Figure [l5j combining the C-AGB and x-AGB photo¬ 
metric classes into the C-rich category. In both galaxies, 
the number of a-AGB stars decreases from Sequence 2 to 
1, and the number of C-rich stars increases in the same 
luminosity range, hinting that these new C-AGB stars 


originated from the Sequence 2 a-AGB sample. More¬ 
over, the change in the ratio of a-AGB stars to C-rich 
stars from Sequences 2 to 1 is larger in the SMC by a fac¬ 
tor of 1.4, indicating a higher transition efficiency. This is 
consistent with stellar evolution models, which predict a 
higher dredge-up efficiency at low metallicities due both 
to less available free oxygen and because the depth of the 

dredge-up events increases at low metallicities. _ 

Stellar evolution mode ls (e.g., Karakas et al.| [2002 
Marigo et al. 2008, 2013) also predict that the dredge- 


up efficiency rapidly decreases at lower masses. The late 
transition from O- to C-rich among the a-AGB stars sup¬ 
ports this prediction, suggesting that they require more 
dredge-up events to reach C/O > 1 compa red to their 
more massive C-AGB counterparts. Karakas (2014) com¬ 
pute that C/O in an LMC-like 1.5 M Q star (Z = 0.007) 
will increase by ~0.2-0.25 after each thermal pulse. With 
an initial C/O near 0.4, LMC a-AGB stars thus require 
a few thermal pulses to become C stars. This is reflected 
in the low fraction (23%) of C-rich a-AGB stars in the 
LMC. 


According to the same Karakas] ([2014]) models, an 

1.5 Mq), becomes C-rich 


SMC-like star (Z = 0.004, M = 
after a single pulse, with an increase in C/O of 1.5. We 
would therefore expect about half of the a-AGB stars to 
be C-rich if they are all at or above the O-to-C transition 
mass, and indeed this is exactly what we find. This in¬ 
dicates that there are few stars in the SMC with masses 
below the O-to-C transition that have reached the final 
stage of the AGB phase. 

4.4. Long Secondary Periods 

Stars with long secondary periods (LSPs) lie on pulsa¬ 
tion Sequence D, the cause of which is unknown. With 
the collection of data we have here, we have found some 
clues to the origin of the LSP: (1) Sequence D is domi¬ 
nated by O- and a-AGB stars (80%; Table [6]); (2) DPRs 
on Sequence D are similar to those on Sequences 1 and 
2; (3) unlike on other Sequences, the ratio of C/M stars 
on Sequence D is almost the same in both galaxies (0.39 
in the LMC; 0.46 in the SMC), despite the difference in 
metallicity; and (4) there are almost no x-AGB stars on 
Sequence D. 

The lack of x-AGB stars on Sequence D may be par¬ 
tially caused by circumstellar extinction; a significant 
fraction are not detected as variable in O GLE, b ut ar e 
recovered as variable in the infrared ( Riebel et al. 12015 ). 
If this is the case, then any x-AGB stars with LSPs are 
dustier than those found on Sequence 1. Of the a-AGB 
stars that show significant dust production (based on 
whether they are detected at 24 /im), 19% (LMC) and 
32% (SMC) are found on Seq uence D. This is consistent 
with Wood & Nicholls (2009), who find that the LSP is 
a ssociated with dusty mass ejection. 

Wood (2015]) show that the primary periods (i.e., the 


seconddongest period) of most LSP stars lies near Se¬ 
quences 2 and 3, as expected if most of the LSP stars are 
O- or a-AGB. The preference for stars with low initial 
masses and O-rich chemistries on Sequence D is a com¬ 
pelling hint for the nature of the LSP and may point to a 
connection between the LSP and the initial onset of dust 
production. A higher fraction of the a-AGB stars reside 
on Sequence D in the SMC, suggesting that metallicity 
may also play a role in the origin of the LSP. 
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Fig. 16. — LMC Spectral classifications compared to the J — [8.0] 
color. The upper and lower panels show the spectrally-classified C- 
rich and O-rich stars, respectively. Points outlined in black were 
photometrically classified as a-AGB stars. The solid line marks the 
best photometric division between C-rich and O-rich stars. 

4.5. Implications for Photometric Classification 

Even though the a-AGB stars have distinct J — [ 8] col- 
ors frornthe O-AGB and C-AGB populations, the |Ciom 
et al. (2006) J — K$ classification scheme works rela¬ 
tive ly well in the MCs from a statistical perspective (Sec- 
tion|2.3.1 ). However, this scheme fails when the goal is to 


classify the most evolved, dust-producing sources, includ¬ 
ing both the a-AGB stars and the more massive O-AGB 
stars with high DPRs (Fig . [l6). The results for M31 de¬ 
scribed in Boyer et al. (201o|suggest that the fraction of 
misclassified O-rich stars will also increase with metal- 
licity due to strong water absorption in late-type O-rich 
giants. 

The black line in the lower panel of Figure [16] shows 
the best J — [8] separation in the LMC, resulting in the 
correct classification for 93% of the C stars and 92% 
of the O-rich stars. The misclassifications may be par¬ 
tially explained by non-simultaneous observations be¬ 
tween 2MASS and Spitzer , which can cause artificially 
blue or red colors owing to variability. The James Webb 
Space Telescope (JWST) will include filters covering 0.7- 
25.5 /im, allowing for near-simultaneous imaging over a 
wide wavelength range. JWST will also have a filter cen¬ 
tered on the 10-/im silicate feature present in most dusty 
O-rich stars, which may be a better choice than the IRAC 
8.0-/im filter for separating dusty M and C stars. 

Even though the statistical separation of C and M stars 
is only marginally improved by including the 8-/im pho¬ 
tometry, we recommend using the J — [8] color to distin¬ 
guish evolved TP-AGB stars (especially the a-AGB stars) 
from their less evolved counterparts when pulsation in¬ 
formation is not available. This separation is impossible 


in J — Ks and in various combinations with the other 
three IRAC filters. We show the J — [8] color distri¬ 
bution for stars on each pulsation sequence in Figure p%| 
there is overlap depending on the exact definition of each 
pulsation sequence, but it is clear that stars redder than 
the solid line are dominated by evolved AGB stars on se¬ 
quences 1 and 2 (a-AGB, C-AGB, x-AGB, and massive 
O-AGB stars). Without this or a similar color cut, it 
is difficult to identify the highly evolved, low-mass stars 
(the a-AGB stars) without either the pulsation informa¬ 
tion or complete spectral energy distribution fits for es¬ 
timating the DPR. 

The [8] — [24] color (Fig.[l4|) is also a reliable way to sep¬ 
arate the evolved a-AGB stars from less-evolved O-AGB 
stars. However, it results in an incomplete sample owing 
to the limited sensitivity at long wavelengths. JWST will 
be approximately 3-4 times more sensitive than Spitzer 
at 24 /im. 

4.6. The Cause of the Red J — [8] colors 

Bll first noticed the a-AGB stars by their red J — [8] 
colors. Their colors must be caused by the flux in the 
8-/im filter since the a-AGB stars do not stand out in 
other color combinations involving J. There are two pos¬ 
sible spectral features that might affect this wavelength 
range: SiO absorption at 8 /im and silicate dust emission 
at 10 /im. The a-AGB stars have modestly high DPRs, 
confirming the presence of dust. Dust can cause both 
strong silicate emission and/or veiling of the SiO molec¬ 
ular band, and both scenarios would increase the J — [8] 
color. Ten of the LMC a-AGB stars have been observed 
with the IRS and their AGB nature confirmed. Of these, 


nine show strong si licate emissio n (Fig . 18). 

In dustless stars, Sloan et al. (2008) sEowed that the 
strength of SiO absorption decreases with metallicity. If 
SiO is the dominant factor causing the red J — [8] color, 
we would therefore expect redder J — [8] colors for the 
SMC a-AGB stars. Figure J7] shows the distribution in 
a-AGB colors, and the SMCJcolors are, in fact, slightly 
bluer , on average. We therefore conclude that the silicate 
feature is the dominant cause of the red J — [8] color for 
O-rich a-AGB stars. The J — [8] colors of the massive 
dusty O-AGB stars are similar to the a-AGB colors (see, 
for example, Fig. 12). It is therefore reasonable that the 
red colors are caused by the same mechanism, i.e., the 
silicate feature. 

Bll speculated that the red J — [8] colors of a-AGB 
stars might be caused by an unusually shaped silicate 
feature. The silicate features of these nine examples are 
not remarkable compared to other O-rich AGB stars in 
the S AGE-Spec program (|Woods et al.|2011||Ruffle et al. 
2015), so we rule out an unusual feature shape as a cause 
of the red colors of a-AGB stars. Instead, the red colors 
are more likely caused by the overlap between the wing 
of the silicate feature and the 8-/im filter transmission 
(Fig.®. 

Dusty C-rich a-AGB stars would also be red in J — 
[8] because of the contribution of the dust continuum 
of amorphous carbon dust at 8 /im; dustier stars would 
exhibit redder colors. Compared to other C-AGB stars, 
the C-rich a-AGB stars have low DPRs, so their colors 
at the blue end of the C-AGB branch in Fig. [I] is to be 
expected. 

C-rich and O-rich a-AGB stars are similarly indistin- 
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J-[8] J-[8] J-[8] J-[8] J-[8] 


Fig. 17. — Distribution in J — [8] color for stars on each pulsation sequence. The most evolved stars, particularly the low-mass (a-AGB) 
stars, can be identified as those redder than the solid line that intersects the upper and lower axes at J — [8] = 1.3 mag and J — [8] =3 mag. 



Fig . 18.— IRS spectra of a-AGB stars in the LMC (Kemper et al. 
2010), in normalized flux density. Nine show silicate emission at 
10 jum. The 8-/xm IRAC filter (grey line) clips the wing of the 
silicate feature, causing red J — [8] colors. The stars with the 
strongest silicate features (in the lower three panels on the right) 
show J — [8] > 2 mag. 


rich, we now refer to the Bll aO-stars as simply a-AGB 
stars, and find that they are evolved, low-mass dusty 
AGB stars. We reach this conclusion by interpreting the 
following observations: 

• Their optical spectra show a high fraction of both 
C- and O-rich chemistries despite having J — K$ 
colors that are typically assigned to O-rich stars. 
There is a higher fraction of C-rich a-AGB stars in 
the more metal-poor SMC (50%, compared to 23% 
in the LMC). 

• Light curves indicate they are pulsating primarily 
in the fundamental mode and first overtone, sug¬ 
gesting they are at the very end of their evolution. 
Among those pulsating in the fundamental mode, 
the a-AGB median pulsation masses are lower than 
those of the O-AGB and C-AGB populations. 

• We identified (likely intrinsic) S-type stars in the 
SMC, and find that they comprise 9% of the a-AGB 
population. These pulse primarily in the first over¬ 
tone and cluster near [3.6] = 11.3-11.8 mag, sug¬ 
gesting that stars with similar properties will tran¬ 
sition to C-rich as they move to the fundamental 
mode. 

• Their dust-production rates are consistent with 
low-mass O-AGB stars. 


guishable in color combinations that include the AKARI 
ll-/rm flux since both C-rich and O-rich dusty stars have 
strong excess at 11 fi m. 


5. CONCLUSIONS 


Using Spitzer SAGE data, Boyer et al. (2011]) noted a 
group of thermally-pulsing AGB stars in the Magellanic 
Clouds which classify as O-rich based on their J — K$ 
colors, but that show redder J — [8] colors than the bulk 
of their parent population. That work dubbed these the 
anomalous O-rich (aO-) AGB stars. We have gathered ev¬ 
idence from multiple perspectives to ascertain the nature 
of these stars and conclude that they are low-mass AGB 
stars that are straddling the mass limit where dredge up 
is sufficient to transform them into carbon stars. Since 
a high fraction of the aO-AGB stars turned out to be C- 


• Their colors match stellar evolution tracks for low- 
mass AGB stars. 

• Their metallicities are low compared to the O-AGB 
stars, implying that they originated from an older, 
less chemically-evolved, stellar population. 

The results here suggest that a-AGB stars can be easily 
selected photometrically using the J — [8] color, which 
is more reliable than J — K$ for the most evolved stars, 
and particularly for the dust-enshrouded O-rich stars. In 
addition, surveys searching for S-type stars may be wise 
to search among the brightest a-AGB stars. 

With the compiled rich data set for the Magellanic 
Clouds, we have noted some general properties of the 
a-AGB stars in addition to the evidence listed above. 
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First, their dust-production rates are not strongly depen¬ 
dent on their periods, at least compared to other O- and 
C-AGB stars pulsating in the fundamental mode (Mira 
variables). Second, they are over represented among 
stars with a long secondary period. Third, their ef¬ 
fective temperatures are similar to high-luminosity O- 
rich evolved stars. Finally, even those that will even¬ 
tually transition into C-AGB stars will produce signif¬ 
icant quantities of O-rich dust prior to this transition. 
More massive stars will become C-rich before the onset 
of strong dust production (except the most massive hot- 
bottom-burning AGB stars, which will remain O-rich). 

In the LMC, models suggest the a-AGB stars have ini¬ 
tial masses up to 1.5 M©, and up to 1.25 M e in the 
SMC. The transition mass would therefore be a bit be¬ 
low this limit. In both of the Magellanic Clouds, there 
is a peak in the star-formation history th at corresp onds 
to initial stellar masses of 1—1.3 M 0 (|Harris & Zaritsky 
2004 2009). The evidence presented here suggests that 
this star formation event is the origin of the a-AGB stars, 
and may have been caused by an encounter between the 
LMC and SMC «2.5 Gyr ago. 
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